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Abstract 

The technique of Caffarehi and Silvestre |T], characterizing the fractional Laplacian 
as the Dirichlet-to-Neumann map for a function U satisfying an elliptic equation in the 
upper half space with one extra spatial dimension, is shown to hold for general positive, 
non- integer orders of the fractional Laplace operator, by showing an equivalence between 
the norm on the boundary and a suitable higher-order seminorm of U . 

1 Introduction 

In [l], Caffarelli and Silvestre showed that, for < 7 < 1, the fractional Laplacian (— A)"^ 
of a function / living on can be understood as the Dirichlet-to-Neumann map for a 
function U living on the upper half-space , where U coincided with / on M", and 
U satisfied a particular 2nd-order ehiptic equation. The function U was thus a suitable 
extension of /. The extension was proved in two separate ways: firstly through an analysis 
of the Poisson kernel, and secondly through an energy equality where the H'^ seminorm 
of / was shown to be equivalent to a suitable Dirichlet-like energy for U, taking the form 
y^-'^^\VU\'^dxdy. 

In this note, we generalize the energy equality of Caffarelli and Silvestre to show that the 
fractional Laplacian of any positive, non-integer order can be represented as a higher-order 
Neumann derivative of an extended function U, where U satisfies a higher-order elliptic 
equation. 

To illustrate the technique, we first show that in the case 1 < 7 < 2, the fractional Lapla- 
cian (—A)''' can still be represented as a suitable Neumann derivative for the solution of a 
higher order equation, and subsequently we generalize this to all positive, non-integer values 
of 7. Chang and Gonzalez [2] showed that the extension has an interesting interpretation 
in terms of scattering theory, and, in particular that the following equation 

AU + -Uy = (1) 
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(where a = 1 — 27) holds for the extended function [/, and still holds for higher orders, 
so long as 7 < ^ (for n odd, and all higher orders for n even). The equation ([1]) plays an 
interesting role in the extension; specifically, we will see that it holds even when 7 > ^, and 
indeed does so for all noninteger 7. 

We conclude the paper with a new proof for strong unique continuation for fractional- 
harmonic functions of higher non-integer order. The proof uses both our extension technique 
to characterize fractional harmonic functions as boundary data for a very high order elliptic 
equation. We then use a technique by Garofalo and Lin 0], which uses a variation on the 
Almgren frequency formula to demonstrate unique continuation. 

The extension technique of Caffarelli and Silvestre has previously been adapted to an- 
alyze fractional-order powers of Schrodinger operators by Stinga and Torrea [SJ, although 
the overall order of the operator remained 7 < 1. Fall and Felli |3] have recently used 
the Almgren frequency formula to analyze unique continuation for fractional-order elliptic 
equations, although their analysis is also restricted to 7 < 1. 

2 The model case: 7 = | 
2.1 

First we discuss the extension in a special case, which illustrates the main point of the 
argument without the complexity of notation we need for more general cases. In what 
follows, 7 = |, and a = 1 — 27 = —2. 

Theorem 2.1. For functions U G Ty^'^(M"'^^) satisfying the equation 

A^U{x,y) = (2) 

on the upper half space for {x,y) G M" x M_|_, where y is the special direction, and the 
boundary conditions 

U{x,0) = fix) 
Uy{x,0) = 

along {y = 0} where u{x) is some function defined on H"'{M.'^) we have the result that 

i-A)lf{x) = Cn,^-^AUix,0) (3) 

Specifically, 

[ \^?\f{0?d^ = Cna I \AU{x,y)\^dxdy (4) 
Proof. Taking the Fourier transform in the x variable only on the energy term Au we get 

-\(\'u{^,y) + Uyy{C,y) 

And so minimizing the energy corresponds to minimizing the integral 

/ / \-\^\^uity) + Uyy{C,y)\^dCdy. 
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Integrating by parts, we see that for each value of ^, U solves the ODE 
Let 6 G VF^'^fM , ) be the minimizer of the functional 



J{4>)= f {f{y)-m?dy 



among functions satisfying the conditions <j){Q) = 1, i?!>'(0) = 0. Thus cf) solves the ODE 

<t>{y)-2cP"{y) + <f>""{y)=0 

with appropriate boundary conditions and we see that U{^,y) = /(O^^dClz/) is a good 
representation for U. 

Calculating, we see that 



/ 



{Mjfdxdy = Cn I \-\i\^U +Uyy^ 



d^dy 

Cnj \-\efmmy) + \i?fm"myi\d^dy 
Cnj \^\'\fm\-m+<t>"{y)f^^di 



and hence the energies are identical up to a constant. 

The Euler-Lagrange equation for the left hand side above is simply the bi-Laplace equa- 
tion, while for the right hand side it is the fractional harmonic equation of order 7, and the 
rest of the result follows. □ 

2.2 

We remark that, using the equation 

y 

and the boundary conditions, we can perform a brief calculation 
/ {AUfdxdy = - [ {AU)^dx - [ V{AU)-VUdxdy 

Jy>£ Jy=e Jy>e 

= [ (AU)^dx-2 [ v(^YvUdxdy 
Jy=s dy Jy>^ \y J 

Jy>£ r y 

Jy=e dy Jy=e 2 2y Jy^, 



y>e y'^ y"^ 



dxdy 
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We take the limit as e — )• and use that Uy{x, 0) = 0, to recover 



If {AUfdxdy = lim(- [ \V^f?dx-f ^^^dxdy] 

for solutions U of our extension problem. This relationship gives a relationship between 
the new energy for the extension in this case, and the straightforward extension of the 
old energy, J y^~^'''|VC/|^, which is infinite in this case. If we take the point of view of 
considering the extension problem as a version of the scattering problem a la ]2], we have 
recovered our energy, J {AU)'^dxdy, as the equivalent of the finite part, of the energy integral 
J y~^|VC/p (see, e.g., O Proposition 3.2]). A similar calculation is possible in the general 
of non-integer 7, yielding many more boundary terms which blow up to infinity at different 
rates. 



3 1 < 7 < 2 

In this case, the argument is precisely analogous to the previous section, except that, like 
Caffarelli and Silvestre, we shall use a weighted seminorm. To be precise, we attach the 
weighted measure y^dydx to our Sobolev spaces, and consider energy minimizers with re- 
spect to this measure on the upper half space of an appropriate energy. Here, we take 
6 = 3 - 27. 

To construct the appropriate energy in this space, we introduce the following operator, 
which is a variant of the Laplacian adapted to the measure, whose virtue is that in the 
weighted space it behaves under integration by parts just as the regular Laplacian does in 
an unweighted space. 

AbU = AU + -Uy 

y 

gives us the desired relationship: 

[ (V$ • V'^)y^dydx = - [ ^ lim ( y^^] dx - [ ^{Ab^)y^dydx 

Clearly, the appropriate 2nd order seminorm for our space is 



/ 



y^\AbU\^dydx 

ri + l 



Our space will be equipped with the norm 

2 — A ii2 — \ / ii2 — ii2 



Our main result is now that 
Theorem 3.1. For functions U G W"^'"^ {MJ^^ , y'') satisfying the equation 

AlU{x,y) = Q (5) 
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on the upper half space for {x,y) G M" x where y is the special direction, and the 
boundary conditions 



U{x,0) 
limy^C/j,(x,0) 







(6) 



along {y = 0} where f{x) is some function defined on H'^{W^) we have the result that 

{-Ay fix) = Cn,^ lim/|-Abt/(x,y) 
y^o ay 

Specifically, 

\i?''\i\0?di = Cn,J' ^y'>\A,U{x,y)\^dxdy (7) 



/ 



n+l 



Proof. Existence and uniqueness of a solution is guaranteed by the usual considerations. 
Taking the Fourier transform in the x variable only on the equation A^u = we get 



\itU - (^ICP + '-^)Uy + (-2|eP + '-^)Uyy + -Uyyy + ^ 



yy ^ y ^yyy ^ ^yyyy 







which is a 4th order ODE in y for each value of ^. Let (f) G PF^'^(R+,y^) be the minimizer 
of the functional 



m= I y\^"{y) + -<t>'{y)-m?dy 
Jr+ y 



among functions satisfying the conditions 0(0) = 1, 0'(O) = 0. Thus cf) solves the ODE 

y y6 yJ. y 

with appropriate boundary conditions and we see that U{^,y) = f{^)4'{\^\y) is a good 
representation for U. 

Calculating, we see that 



/ 



y\Ai,Ufdxdy 



Cn 



-\^\^U+-Uy + Uyy 



d^y'dy 



= CnJ \^\V{Ofi-m-^^^'iy) + ^'\y)f0Td^ 

= c^m J \in\o\'d^ 

and hence the energies are identical up to a constant. 

The Euler-Lagrange equation for the left hand side above is simply the modified bi- 
Laplace equation, while for the right hand side it is the fractional harmonic equation of 
order 7, and the rest of the result follows. 

□ 
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Remark: It is not surprising (and easily verified) that solutions to AU + jjUy = 0, with 
appropriate Dirichlet conditions, satisfy the equation A^f/ = 0. From the scattering theory 
(see, e.g., [2]), we know that limy^Qy^Uy = for this equation, and thus from the unique- 
ness of solutions to the modified biharmonic equation we can see that the function of the 
extension described here, and the one in [2], are identical. 



4 The General Case 



4.1 



The general case follows on a similar theme, taking progressively higher powers of the 
weighted Laplacian A;,. In setting our boundary conditions, we take our cue from [2], 
whence we learn that, when 7 < § the extension function satisfies 

1-27 
AU + ^-Uy = 

y 

and furthermore that, ifm<7<m + l, U has a series expansion in y that has only even 
integer powers, up to the power y'^^. Being inspired by the belief that these considerations 
should hold for all non-integer values of 7, we set boundary conditions for the y-derivatives 
of U in the following result. 

Theorem 4.1. Let 7 > 6e some non-integer, positive power of the Laplacian. Let m < 
7 < m -h 1, orm= [7], and b{-f) = 2m + 1 - 27. For functions U € W"'+'^''^{M.l~^^ ,y^) 
satisfying the equation 

A^+^U{x,y)=0 (8) 

on the upper half space for {x,y) € M" x R_|_, where y is the special direction, and the 
boundary conditions are that U{x,0) = u{x) along {y = 0}, and, furthermore, that for 
every positive odd integer 2k + 1 < m + 1, we have limj^^o y^ ^Qy2k+i (^1 0) = 0; where f{x) 
is some function defined on H'^{W^). For even integers, we specify the relationship 

-(x,0) = (A^C/(x,0))n 



Qy2k^ ' ^ ^ ^-^■-'"^^JLi27-4(j-l 
Then we have the result that 

d_ 

y->o " dy 

Specifically, if m is odd. 



i- Ay fix) = Cn,-y lhny—A^U{x,y) (9) 



m+l 



\^\^^\KO\^dC = Cn,j y'lA,^ U{x,y)\^dxdy (10) 



n + l 
+ 



and if m is even, 



r I m-f-l I 

!el'^l/(Ol'de = C„,^ / y''\VA[—^U{x,y)\^dxdy (11) 
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Proof. The proof follows along the same lines as before, although some of the auxiliary 
calculations can be a bit complicated, and are safely left to a lemma that we prove later. 
As in the previous cases, we consider a representation for the solution given by U{£,,y) = 
If m is an odd integer, then we let (/> be the minimizer of 

subject to (/)(2fc+i)(o) = for ah odd 2k + I < 2m + 1, and 0(0) = 1, with (/>(2fc)(0) = 
Y['j=i 2j-l{j-i) • ^ even integer, then we let 



W[(l),y] ( yQy^ g2y^ 



and take 

J(0) 



({W[^,y])^ + (^^W[(P,y]^ jdy 



The equivalence of the energies then follows precisely as in the argument in Theorem 13. 1[ 
Only the calculation of the Euler-Lagrange equation for our energy functional remains. 
It is readily seen that A™+^C/ = 0, but it is necessary to show that boundary terms of the 

form 4„ limy^oy''^^^^^p^A'^-''U{x,y)dx = 0. In fact, = o. The calculation 

is technical, and we leave it to Lemma l4.4i □ 

Lemma 4.2. Functions U satisfying the equation 

1-27 

AU+ -Uy = 

y 

Ar+^c/ = 



also satisfy 
and W = A^U satisfies 



y 



for any k < m + 1. Here m, b are defined in terms of j as in the statement of Theorem\4.1 



Proof. The proof shall be an induction on m. When m = 0, this is a restatement of the 
definition of Af,. Notationally, let us write a = 1 — 2^. 

Suppose now, for general m, that AU + fUy = 0. Let W = A^f/. Then 



AW = ^[^U, 

and by the inductive hypothesis the result follows. □ 
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Remark It is no surprise that, having carefully set the boundary conditions to coincide 
with the function U from scattering theory [2], that our energy minimizer, satisfying the 
same equation as the U of the scattering theory, would be exactly the same function by the 
uniqueness of energy minimizer s. 

We shall now need a partial converse to the last lemma. In the case of 7 < ^, it is 
already known, and a consequence of the scattering result. 

Lemma 4.3. Solutions of the equation A^^^U = with boundary conditions given in 
Theorem \4-l\ also satisfy the equation 



1 - 2t 

AU+ -Uy = (12) 

y 

Proof. We may assume 7 > 1, because otherwise our statement is that of the original 
Caffarelli-Silvestre extension. 

In any subdomain strictly bounded away from the line y = 0, the equation (|12p is a 
classical 2nd order equation with bounded and regular coefficients, and hence we can solve 
the equation with boundary conditions given by the energy minimizer U; let's call this new 
solution U'. A result of Lin and Wang [6] gives C^'" up to the boundary regularity for 
solutions of the equation (I12p . Hence we may approach by taking domains that approach 
ever more closely to the boundary, reapplying our argument, and passing to the limit. In 
the limit, it is no surprise that the limiting value of U' is the same as the original U, since 
the two both satisfy {Ah)"^+^U = with the same boundary conditions. □ 

We now have the tools to complete the proof: 

Lemma 4.4. For minimizers U of our energy satisfying the preconditions of Theorem \4-l[ 
we have that 

a(A^^(x,0)) _^ 
dy 

for k < m. 

Proof. Let W = A^C/. Then W solves A^~^^~''W = 0, and hence satisfies the equation 

a + 2k^ 

y 

By the maximum principle, W{x, 0) + Cy^ is locally a supersolution for C > sufficiently 
large, and W{x, 0) — Cy'^ is a subsolution, again for C > sufficiently large. Hence, we can 
conclude that ^(x, 0) = 0. □ 

As a consequence of our results, we have the following trace inequality: 

Corollary 4.5. For U £ ^^"+l(y^ satisfying the given boundary conditions, we have 

that the trace of U on the boundary, which we call f , satisfies 

H7(]Rn) < C.y\\y^ D^^^U\\^2l^^^+i) 



AW + — Wy = 



where is either A, ^ , or VA,^ , whichever is appropriate. 



m+l 

The proof follows from the energy argument of Theorem 14. 1^ since that theorem uses 
energy minimizers. 
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5 Strong unique continuation 



We show an analogue of Almgren's frequency formula holds for the fractional Laplacian, 
which allows us to deduce strong unique continuation in the style of Garofalo and Lin [3]. 

We begin by showing that such a formula holds for balls centered on the interior of a 
domain where the equation A™^^C/ = holds. The calculations along the boundary are not 
so very different, and the boundary terms are absorbed Let Uk = {Ah)^U, and we consider 
the behavior of certain integrals of U with respect to balls centered at a particular point. 
Let 



and 



Then we let 



III n 

D{r) = y' d^^'^l' + UkUk+i) dxdy 

m „ 

H{r) = / y'iUkfdS. 



and we get the following result 

Theorem 5.1. There exists a number A{m,b,n) > such that if (Af,)"^~^^U = 0, then 
e^^N{r) is uniformly bounded for r < 1. 

In the proof, an important calculational tool is the following modification of an identity 
originally due to Rellich [7], which appeared in [T]. 

Lemma 5.2. Suppose AbW{x,y) = V{x,y). Then 
rj y^{\VW\^ -2W^)dS = {n-l + b) I y^VWfdxdy-lf y\X-VW)Vdxdy 

'J OBr ^ Br ^ Br 

Proof of Lemma We apply the divergence theorem to the following term. 

V • (y'^^X - y\X ■ VW)VW^ 

The boundary term is precisely the desired left hand side, and after some calculation we 
find that the right hand side matches up with the given term as well. □ 

Another extremely useful estimate is the following calculation, which allows us to bound 
the behavior of an integral in the interior of the ball in terms of its behavior on the boundary. 

Lemma 5.3. Suppose AijW{x^y) = V{x,y). Then 

[ y^W^dxdy <c(r [ y^W^dS + / V^dxdy) 

J Br V JdBr JBr J 
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Proof of LemmalEE Notice that AbW^ = 2\VW\'^ + 2WAbW = 2|VPFp + 2WV. We 
consider the integral 

/ y^AbiW"^) (r^ - dxdy = j y^V{W^) ■ i^ir^ - \X\'^)dS - j y^2WS/W ■ {-2X)dxdy 

J Br OBr Br 



4 / y^WiVW ■ X)dxdy 

J Br 



= 4/ yVrdS- 4WV ■ (y''WX)dxdy 

JdBr J Br ^ ^ 

= Ar I yVdS - 4 / y^W [(n + 1 + h)W + {VW ■ X)] dxdy 

J OBr ^ Br 

Taking the last line of this calculation and comparing it with the second, we see that 

/ y^W{VW ■ X)dxdy = l r [ y^W^dS -(n + l + b) I y^W'^dxdy 
and hence that 

/ y^\VW\^dxdy + {n + l + h) I y'^W'^dxdy = r f y^W^dS - f y^WVdxdy 

J Br ^ Br ^ OBr ^ Br 

Applying Cauchy's inequality, we get the desired result. □ 

With these two estimates in hand, we can proceed to the proof. The general idea of the 
proof is the following: the map r i— N{r) is continuous in r, and so we can divide it into 
the set where N{r) is large, and where it is small. We are only interested in proving our 
estimate when N{r) is large, for the rest of the time it is already known to be small. 

Proof of Theorem \5.1l We confine ourselves to the study of the set where N{r) > 1. 
We consider the expression logA^(r), and take its derivative, whence we get 

1 _^ D'{r) H'{r) 



r D{r) H{r) 
Calculating and using the lemmata above, we find that 



D'{r) 



E 



/ ((vc/fc)2 + UkUk+i) dS 

JdBr 
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We examine first the gradient square term 

n-l + b 



dBr 



y\VUkYdS 



B, 



y\VUkYdxdy + 2 y" 

I dBr 



r J Br 

n-l + h 



y\X ■ VUk)Uk+idxdy 



B, 



y\VUkYdxdy + 2 y 

J dBr 

n+l + h 



dr 



dS 



dr 



dS 



2/ y'U.U.^.dS + l'^^^ I 

JdBr JBr 



y UkUk+idxdy 



2 

+- 

r 



y'Uk+iX ■ VUkdxdy 



Br 



Hence if we write 
then we have 



Dk{r) = I y'(|VC/fc|2 + C/fc?7fc+i) 



dydx 



n-l + h 



Dk(r) + 2 



dBr 



dUk 



dr 



dS-2 y'UkUk+idS 

I dBr 



+- 



n + 3 + b 



1 y^UkUk+idxdy + - j y'^Uk+iiX ■ VUk)dxdy 

J B-p ^ J Br 



We now show that the last 3 terms of this expression are bounded in terms of H{r) 
The first term is a straightforward appHcation of elementary inequalities: 



dBr 



y'UkUk+idS 



< 



dBr 



y' W + \Uk+i\'') dS 



When we sum over all values of k, this is clearly bounded by 2H[r). 

The second term is bounded in a similar fashion. First, an elementary inequality: 



y UkUk+idxdy 



Br 



< 



followed by iterated use of Lemma [ 

- [ y'Uidxdy <c( [ y'UidS + / y'U^^.dxdy) <■■■< CH{r) 

^ JBr \JdBr JBr J 

The third term we control in terms of Dk{r) and H{r). 

-j y'Uk+i{X-VUk)dxdy< - f y'ul^^dxdy + j y'\VUk\^dxdy 

^ t/ B-p ^ J Br J Br 

Then by iterated use of Lemma [5T3l we can bound ^ y^U^_^^dxdy by H{r), while the sum 
over k of y^\VUk\'^dxdy is clearly bounded by D{r) + H{ 



r . 
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Hence, we conclude that 



D'{r) _{n-l + b) Efc IdBr y' 



dUk 



dr 



D r) 



D{r) 



dS 1 

— + ^^{0{H{r)) + 0{D{r))) 



Now, 0{D{r)) / D{r) is bounded, and in the range we consider, where r < 1 and N{r) > 1, 



so IS 



H(r) 



We now turn our examination to the study of -jj^ ■ A straightforward calculation yields 



that 



dr 



We note that Dk{r) = Jq^ y^Uk^^dS through a simple integration by parts. 
Hence, in the range where r < 1, N{r) > 1, we can estimate that 



log(iV(r))' 



1 ^ D'{r) H'{r 



r D{r) H{r) 



0(l) + 2 



D{r) 



V 



D{r) 



The second term is non-negative by Holder's inequality, and hence 

log(iV(r))' > -C 
for some suitable constant C, which establishes the theorem. 



□ 



Every step of the proofs contained above also apply along the (suitably smooth) bound- 
ary of a domain where U solves A™^^[/ = 0, provided the all the appropriate normal 
derivatives are 0, and hence, if we define, about points (x,0) G M""*""^, 

m „ 

D{r) = V / (|VC/fc|2 + UkUk+i) dxdy 



and 



. and 



m „ 

H{r) = ^ / 

N{r) = r 



{y>0} 

D{r) 
H{r) 



y\UkfdS 



we get the following analogous result: 

Theorem 5.4. The result of Theorem \5.1\ still holds if we consider functions U which solve 



Af'+^C/ = m the upper half space with S-A'^U = OforO<k<m. 



12 



Following the proof of Garofalo and Lin [5] (pp 256-257), that such an integral is 
bounded, along with our extension result, suffices to provide a proof of strong unique con- 
tinuation, that is. 

Corollary 5.5. Let f(x) be a function on M" such that 



in some smooth domain Q C M", with f{x) defined suitably on \ fl, so that f G i7'^(]R"). 
Then if f vanishes of infinite order at some point xq S f^, that is to say for sufficiently 
small r, we have 



for every A; G N, then f is uniformly in Q. 

Funding 

This work was supported by a NSF Mathematical Sciences Postdoctoral Research Fellowship 
(DMS Award number 1103786). 

Acknowledgements 

1 am indebted to several people for their encouragement and advice in the preparation of this 
paper. Fang Hua Lin suggested that a suitable variation of Almgren's frequency formula 
would hold for this extension and suggested viewing it as an elliptic system, leading to the 
discussion in ^ Qing Jie suggested the characterization of our new energy as the finite 
part of the divergent energy integral J y"'\VU\'^dxdy, leading to the discussion in §2.21 Alice 
Chang asked the question that led to this work, and provided much moral support. 

References 

[1] L.A. Caffarelli and L. Silvestre, An extension problem related to the fractional Laplacian, 
Comm. Partial Differential Equations 32 (2007), no. 7-9, 1245-1260. 

[2] Sun- Yung Alice Chang and Maria del Mar Gonzalez, Fractional Laplacian in conformal 
geometry, Adv. Math. 226 (2011), no. 2, 1410-1432. MR 2737789 (2012k:58057) 

[3] Mouhamed Moustapha Fall and Veronica Felli, Unique continuation property and lo- 
cal asymptotics of solutions to fractional elliptic equations. Preprint, arXiv:1301.5119 



[4] Nicola Garofalo and Fang- Hua Lin, Monotonicity properties of variational integrals, Ap 
weights and unique continuation, Indiana Univ. Math. J. 35 (1986), no. 2, 245-268. MR 



833393 (88b:35059) 

[5] C. Robin Graham and Maciej Zworski, Scattering matrix in conformal geometry, Invent. 
Math. 152 (2003), no. 1, 89-118. MR 1965361 (2004c:58064) 



(-A)V = 




(2013). 



13 



[6] Fang Hua Lin and Lihe Wang, A class of fully nonlinear elliptic equations with singular- 
ity at the boundary, J. Geom. Anal. 8 (1998), no. 4, 583-598. MR 1724206 (2001a:35068) 

[7] Franz RcUich, Darstellung der Eigenwerte von Au + An = durch ein Randintegral, 
Math. Z. 46 (1940), 635-636. MR 0002456 (2,56d) 

[8] Pablo Raul Stinga and Jose Luis Torrea, Extension problem and Harnack's inequality 
for some fractional operators, Comm. Partial Differential Equations 35 (2010), no. 11, 
2092-2122. MR 2754080 (2012c:35456) 



14 



